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ABSTRACT 

The nucleus of M82 has been mapped in several 3mm and 1mm lines of CN, HCN, C2H, 0C3H2, 
CH3C2H, HC3N and HOC + using the IRAM 30m telescope. These species have been purposely 
selected as good tracers of photon-dominated chemistry. We have measured [CN]/[HCN]~5 in the 
inner 650 pc galaxy disk. Furthermore, we have detected the HOC+ 1— >0 line with an intensity 
similar to that of the H 13 CO+ 1^0 line. This implies a [HCO+]/[HOC+] ratio of - 40. 
These results corroborate the existence of a giant photo-dissociation region (PDR) in 
the nucleus of M82. In fact, the low [HCO + ]/[HOC + ] ratio can only be explained if the nucleus of 
M82 is formed by small (r< 0.02-0.2 pc) and dense (n~ a few 10 4 -10 5 cm -3 ) clouds immersed in an 
intense UV field (Go^ 10 4 in units of the Habing field). The detection of the hydrocarbons C-C3H2 
and CH3C2H in the nucleus of M82 suggests that a complex carbon chemistry is developing in this 
giant PDR. 

Subject headings: galaxies: individual (M82) — galaxies: nuclei — galaxies: starburst — ISM: 
molecules — ISM: abundances - radio lines: galaxies 
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1. INTRODUCTION 

M82 is one of the nearest and brightest starburst 
galaxies. Located at a distance of 3.9 Mpc, and with 
a luminosity of 3.7xl0 10 L©, it has been extensively 
studied in many molecules. Compared to other pro- 
totypical nearby starburst galaxies like NGC 253 and 
IC 342, M82 presents systematically low abundances 
of th e molecules NH3, CH3OH, CH 3 CN, HNCO and 
SiO l)Takano et alJ 12003^) . Different explanations have 
been proposed to account for this peculiar chemistry. 
Since all these molecules a re related to dust grain chem- 
istry, iTakano et alJ l|2003|) proposed that the formation 
of molecules on dust and/or evaporation to the gas phase 
is not efficient in M82. 

On the other hand, several studies have revealed 
that the starburst has heavily influenced the interstel- 
lar medium in M82 by producing high cosmic rays and 
UV fluxes. A low density ionized component is filling 
a substantial fraction of the volume in M82 (see e.g 
Seaquist et al. 1996). The molecular gas is embedded 
in this component in the for m of warm (Tfr>5 K) and 
dense (n>10 4 cm" 3 ) clouds ijMao et alJl200C|) . Dense 
PDRs are expected to form in the borde rs of these clouds 
JWolfire et al. 1 1199(1 ILord et alJll996T) . IWolfire et all 
(1990) have modeled the CII, Sill and OI emission and 
estimated a UV field of Go=10 4 in units of the Habing 
field and a density of n~10 5 cm -3 for the atomic com- 
ponent. ISchilke et alJ l)1993|) reported observations of 
the CI ( 3 Pi- 3 P ) line deriving a [CI]/CO column den- 
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sity ratio (~ 0.5) higher than that observed in non- 
starburst galaxies (e.g. [CI]/CO ~ 0.15 in our Galaxy). 
They proposed that the enhanced cosmic ray flux sup- 
plied by supernova remnants could produce the enhanced 
CI emission in M82. Recent results suggest that the 
strong UV flux dominates the chemistry of the molec- 
ular gas in M82. IGarcia-Burillo et alJ lj2002D obtained a 
high-angular-resolution image showing widespread HCO 
emission in this galaxy. The enhanced HCO abundance 
([HCO]/[H 13 CO+]-3.6) measured across the whole M82 
disk was interpreted in terms of a giant PDR of 650 pc 
size. 

In this paper, we present observations of a selected set 
of radicals (CN, C 2 H), reactive ions (CO+, HOC+) and 
small hydrocarbons (C-C3H2) which are excellent probes 
of the atomic to molecular transition in PDRs. In partic- 
ular, the [CN]/[HCN] ratio has been successfully used as 
a PDR indicator in regi ons with very different physica l 
conditions in our Galaxy iFuente et alll99.!ll9Mll99fil) . 
The detection of the reactive ion HOC + is almost unam- 
biguously associated to regions with a high i onizing flux, 
eithe r PDRs or XDRs jFuente et all 120031: iRizzo et all 
120031 lUsero et al.M2004j) . Recent works have revealed 
that the abundances of some hydrocarbons, such as c- 
C3H2 and C4H, are an order of magnitude larger in PDRs 
than those predicted b y gas-phase models ijTevssier et alJ 
l27MIPetveraIll2l)0l . 

2. OBSERVATIONS AND ANALYSIS 

The observations were carried out in June and Novem- 
ber 2004 with the IRAM 30 m radiotelescope at Pico 
de Veleta (Spain). We used 2 SIS receivers tuned in 
single-sideband mode in the 1 mm and 3 mm bands. 
The observed transitions are: CN l->0 (113.490 GHz), 
CN 2->l (226.874 GHz), HCN l->0 (88.631 GHz), C 2 H 
l->0 (87.317 and 87.402 GHz), c-C 3 H 2 2 X 2 ^1 1 (85.339 
GHz), C-C3H2 6i, 6 ->5 ,5 (217.822 GHz), CH 3 C 2 H5 fc ->4 fc 
(85.457 GHz), HC 3 N 9^8 (81.881 GHz) and HOC+ l->0 
(89.487 GHz). The line intensities have been scaled 
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to the main beam brightness temperature. The half- 
power beam width (HPBW) of the telescope is 29" at 
85 GHz, 22" at 115 GHz and 12" at 230 GHz. We ob- 
served three positions across the M82 disk: the nucleus 
(RA(2000): 09' l 55 m 51.9 s , Dec(2000): 69°4'47.11") and 
the two peaks in the HCO emission [offsets (+14", +5") 
and (—14",— 5") hereafter referred to as E. and W. knot 
respectively]. The observed spectra are shown in Fig. 1. 
To compare the 1.3mm and 3mm CN and C-C3H2 lines 
and make excitation calculations, we have derived beam 
filling factors from the H 13 CO + interferometric image of 
IGarcfa-Burillo et a.U l|200l . This assumption is justified 
since H 13 CO + is an optically thin tracer of dense gas (n> 
10 4 cm -3 ) and presents a quite uniform abundance in a 
wide range of physical conditions. 

The CH3C2H and C2H column densities have been 
estimated using the LTE approximation. For C2H, we 
have assumed the typic al rotation tempera ture in galac- 
tic PDRs, T rot = 10 K l|Fuente et alJl2003D . In the case 
of CH3C2H, all the lines of the K-ladder are blended. In 
our estimates, we have assumed that 40% of the emis- 
sion comes from the K=0 line and T rot = 20 K. These 
assumptions are based on the CH3C2H observations of 
Sgr B2 a nd Orion which have simila r physical conditions 
to M82 ijChurchwell fc Hollisl Il98^1 . An LVG code has 
been used to derive densities and beam averaged column 
densities for the other species. In these calculations, we 
have assum ed Tfc =60 K inferred from the NH3 lines by 
IWeifi et a.1.1 pOll. 

3. RESULTS 

We have estimated the hydrogen densities in the M82 
disk by fitting with a LVG code the CN and C3H2 
lines. Hydrogen densities between n# 2 ~5xl0 4 cm" 3 and 
n^f 2 ~2xl0 5 cm -3 are derived from both molecules. As- 
suming n# 2 ~lxl0 5 cm -3 , we have calculated the col- 
umn densities averaged in a beam of 29". The CN 
and HCN column densities are quite constant along the 
galaxy disk with values, N(CN)=2.0±0.5x 10 14 cm" 2 
and N(HCN)=4.0±0.5xl0 13 cm- 2 . The [CN]/[HCN] 
ratio is similar to 5 in all positions (see Table 1). 
As discussed in Section 4, this large value of the 
[CN]/[HCN] ratio is only reached in the most heavily 
UV exposed layers of a PDR. The derived C-C3H2 and 
CH3C2H column densities are also quite uniform along 
the disk with values, N(C 3 H 2 )~1.7±0.4x 10 13 cm" 2 
and N(CH 3 C 2 H)~1.0±0.6xl0 14 cm" 2 . These column 
densities a r e in a gree ment with previous estimates by 
lOike et alJ l)2004ft and iMauersbereer et al] (|199H) . Fi- 
nally, we have not detected the HC3N 9— >8 line towards 
any position with an upper limit to the column density 
N(HC 3 N)<2xl0 12 cm" 2 towards the E. knot. As dis- 
cussed in Section 4 the non-detection of HC3N and the 
derived lower limit to the [c-C3H 2 ]/[HC3N] ratio argue 
in favor of a PDR chemistry in M82. 

We have detected the reactive ion HOC + in the three 
selected positions across the M82 disk. Furthermore, 
the intensities of the HOC + 1— >0 lines are similar, even 
larger, than those of the H 13 CO + 1— >0 lines. Assuming 
the same physical conditions for H 13 CO + and HOC + , we 
derive a [HCO+]/[HOC+] ratio of -40 across the 650 pc 
inner disk, which is two orders of magnitude lower than 
that found in galactic giant molecular clouds (GMCs) 
(jApponi k, ZiurvsllT997j) . However, the H 13 CO+ spec- 



tra are derived from interferometric data. In order to 
asses the amount of missed flux (and hence the H 13 CO + 
column densit y), we have compared our res ults to those 
of HCO + from lNguven- Q-Rieu et al.l l)1992[) . which were 
obtained by single-dish observations. From our H 13 CO + 
spectra, we obtain N(H 13 CO+)~5xlO n cm 2 across 
the galaxy. Assuming a 12 C/ 13 C ratio of 89, this 
implies N(HCO+)~4.5x 10 13 cm" 2 . This value is in 
agreeme nt within a factor of 1.5 w ith previous esti- 
mate bv iNguven-Q-Rieu et aD l)1992j) . Thus, even m 
the most conservative case we can conclude that the 
[HCO+]/[HOC+] ratio is < 80 in the M82 disk. Such low 
values of the [HCO + ]/[HOC + ] ratio have on ly been found 
in the galactic reflection nebula NGC 7023 (Fuente et al. 
1 20031) and in the active nucleus of NGC 1068 l|Usero et al. 
2004) and put severe constraints to the ionization degree 
of the molecular gas in M82. 

4. PHOTON-DOMINATED CHEMISTRY IN M82 

iSternberg fc Dalgarnol l) 19951) found from detailed 
modeling that [CN]/[HCN]>1 arises naturally in the sur- 
face layers (A t ,<4 mag) of dense PDRs. Values of the 
[CN]/[HCN] ratio between 1-3 have been found in pro- 
totyp ical dense galactic PDRs ijFuente et al.lll993t 119951 
1996). We have derived [CN]/[HCN] ratios ~5 in all posi- 
tions across the M82 nucleus suggesting that the molecu- 
lar clouds in this galaxy are bathed in an intense UV field. 
To determine the averaged physical conditions of these 
clouds, we have carried out model calculations using the 
plane-parallel PD R model developed b y Le Bourlot and 
collaborators (jLe Bourlot et al71ll993|) . The model in- 
cludes 135 species (HOC+, CH 3 C 2 H and HC 3 N are not 
included) and standard gas-p hase reactions. Adopting 
G ~10 4 ( Wol fire et al. 1119911) and a total hydrogen nu- 
clei density tih = nni + 2 x n# 2 =4x 10 5 cm" 3 , the model 

predicts that [CN]/[HCN] ratios ~5 are only expected in 
regions at A t ,<5~6 mag (see left panels of Fig. 2). This 
implies an important limit to the cloud sizes in the nu- 
cleus of M82. Since the clouds are bathed in a perva- 
sive UV field, the averaged column density of individual 
clouds should be Nh 2 ^W 22 cm" 2 in order to have aver- 
aged [CN]/[HCN] ratios ~5. 

Strong constraints are also derived from the HOC + 
observations. We have measured a [HCO + ]/[HOC + ] ra- 
tio of ~ 40 in the M82 disk. A simple calculation of 
the CO+/HCO+/HOC+ chemical network shows that 
a high ionization degree, X(e~)> 10~ 5 , is required to 
have [HCO+]/[HOC+]<80 llUsero et a,l .120041) . Our PDR 
model shows that this high electron abundance is only 
found at A v <4 mag (see left panels of Fig. 2). This im- 
plies that the averaged column density of the clouds in 
the nucleus of M82 is N(H 2 )<8xl0 21 cm" 2 . Thus, the 
M82 nucleus seems to be formed by small (r~0.02-0.2 
pc) and dense (n~10 4 -10 5 cm" 3 ) cloud s immersed in a 
UV field of G ~10 4 . IMao et all l|2000ft estimated simi- 
lar column densities for the clouds in the inner 400 pc 
disk of M82 by modeling of the 12 CO, 13 CO and C 18 
lines. Low values of the [HCO + ]/[HOC + ] ratio are also 
found in XDRs with high i onization degree, such us ac- 
tive galactic nuclei (AGNs) ljUsero et all2004|) . However, 
the PAH emission is very low in AGNs since the dou- 
bly ionized PAHs produced by X-rays are very unstable 
ijLeach et al.||198 9|). The intense PAH emission observed 
in M82 l)Normand et all 119951 iForster Schreiber et all 
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2003) indicates that X-rays are not the driving agent of 
the molecular gas chemistry in this galaxy. 

5. HYDROCARBON CHEMISTRY 

The hydrocarbon chemistry in PDRs has been a sub- 
ject of increasing interest both from the theoretical and 
the observational point of view. Recent works have re- 
vealed that the C-C3H2 abundance in PDRs is similar 
to that in dark clouds. This is quite surprising taking 
into account that this carbon cycle is easily photodis- 
sociated. In fact, PDR gas-phase models fall short of 
explaining the observed c-C 3 H 2 and C4H abundances by 
an order of mag nitude i|Tevssier et alJl200l iPetv et all 

2004) . This is clearly seen when one compares the [c- 
C 3 H2]/[HC 3 N] ratio in PDRs and dark clouds (see Table 
1). While both species have similar abundances in dark 
clouds, the [c-C 3 H 2 ]/[HC 3 N] ratio is >10 in PDRs. Since 
both molecules are easily destroyed by photodissociation, 
this suggests the existence of an additional c-C 3 H2 for- 
mation mechanism in PDRs. Several authors have pro- 
posed formation process es of small hydrocarbons l inked 
to th e PAHs chemistry ijTevssier et alJl2004t IPetv et alJ 
12004ft . We have detected widespread emission of C2H 
and c-C 3 H2 in the M82 disk. However, we have not de- 
tected HC 3 N towards any position. The derived upper 
limits to the HC 3 N column density in M82 show that the 
[c-C 3 H2]/[HC 3 N] ratio is at least a factor of 10 larger in 
this starburst galaxy than in the prototypical galactic 
dark cloud TMC 1 , and is similar to that found in dense 
galactic PDRs (see Table 1). This result reinforces the 
scenario where the inner 650 pc region of M82 is a giant 
PDR. 

In addition to C2H and c-C 3 H2, we have detected 
CH3C2H in all the positions across the M82 disk. This 
detection is rather puzzling since this molecule is easily 
photodissociated in regions exposed to intense UV fields 
and therefore its abundance is expected to be very 
low in these regions. Furthermore, this molecule is 
not detected in prototypical galactic PDRs, like the 
Horsehead, which is known to be specially rich in c arbon 
compounds l|Tevssier et al.ll2004UPetv et aDl2004ft . One 
possibility is that the CH 3 C 2 H emission arises in a 
population of clouds similar to the galactic GMCs. 
The molecule CH3C2H is quite abundant in GMCs, 
presenting similar abundances to those of chemically 
related species like CH 3 OH and CH 3 CN. The same 
behavior is observed in other prototypical starburst 
galaxies like NGC 253 where [CH 3 C 2 H]/[CH 3 OH]~l 
and [CH 3 C2H]/[CH 3 CN] -8. However, 
rCH 3 C? Hl/fCH 3 OHl>8 and fCH 3 C^Hl/fCH 3 CNl>25 
in M 82 l|Mauersberger et a.1.1 Il99ll iHuettemeister et alJ 
11997ft . i.e., a factor 3-10 larger than in NGC 253. These 
abundance ratios show that the hydrocarbon chemistry 
in M82 is very different from that of GMCs. 

We think that the over-abundance of CH3C2H is 
also related to the enhanced UV flux in this galaxy. 
In fact, CH 3 C 2 H and CH 3 C4H has been detected by 
iCernicharo et all |2001) in the C-rich proto-planetary 
nebula CRL 618. This nebula is specially rich in complex 
carbon compounds, in particular, polyynes (C 2ll +2H 2 ) 
and methylpolyynes (CH 3 C 2n H). They argued that the 
formation of these complex molecules has started by the 
photolysis of acetylene (C 2 H2) and methane (CH4), both 
processes occurring in the inner PDR of the planetary 



nebula. The photo-erosion of PAHs could also contribute 
to enhance the abundance of these molecules by releasing 
small carbon chains and acetylene to the gas phase. 

6. A GLOBAL SCENARIO FOR THE CHEMISTRY IN M82 

The observations presented in this paper show that 
the strong UV flux drives the chemistry of the molec- 
ular gas in the M82 disk. In this Letter, we explore 
whether the PDR chemistry can solely account for all 
the molecular abundances measured in this galaxy so 
far. For this aim, we compare the molecular abun- 
dances in M82 to those of the prototypical starburst 
NGC 253. Molecules such as NH 3 , CH 3 OH, CH 3 CN, SO, 
SiO and HC 3 N are clearly under-abundant in M82 com- 
pared to NGC 253 llWeiB et al.l200l}lHuettemeister et all 
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Garcia-Burillo et alJ 120001 12001ft . All these molecules 



are easily photodissociated and are not expected in a 
PDR. Other species such as CS, CN, HCN and HCO+ 
present fractional abundances ~ a few 10~ 9 in M82 
which are very similar to tho se measured in NGC 253 
ijMauersbereer fc Henkellll989ft . These molecules belong 
to the reduced group of species which have significant 
abundances at A„<5 mag in our PDR model (see left 
panels of Fig. 2). However, our model predicts abun- 
dances a factor of 10 lower than those observed. Cos- 
mic rays and/or shocks could contribute to enhance the 
abundances of these species. 

To investigate the possible effect of the enhanced cos- 
mic ray flux on the chemistry, we have repeated the 
model calculations with an enhanced cosmic rays flux 
of C=4xl0~ 15 s -1 . This is the value estimated by 
iSuchkov et al.1 l)1993ft to account for the physical con- 
ditions of the molecular gas in M82 and is, very likely, 
an upper limit to the actual cosmic ray flux. The en- 
hanced cosmic rays flux does not change significantly the 
value of the [CN] / [HCN] ratio and the ionization degree 
at A,j<5 mag (see right panels of Fig. 2). However, the 
fractional abundances of HCN, CN and CS averaged over 
the region at A„<5 mag are a factor of 10 larger than 
those obtained with the standard value £=5xl0 -17 s _1 . 
Moreover, these values are in reasonable agreement with 
the observational results in M82. Thus, by including the 
enhanced cosmic ray flux we get a better fit of the model 
to the observational results. However, the HCO + abun- 
dance still remains a factor of 10 lower than the observed 
one. 

Some of the underabundant species in M82 like SiO 
are well known shock chemistry tracers. This sug- 
gests that shocks are not the dominant mechanism at 
work in the galaxy disk. The SiO map reported by 
IGartia-Burillo et alJ l|200H) show that shocks are not oc- 
curring in the galaxy disk but in the disk-halo interface 
because of the hot gas coming out of the plane in the 
nucleus of M82. This is consistent with the interpreta- 
tion of M82 as an evolved starburst compared to NGC 
253. While in NGC 253 the molecular chemistry is deter- 
mined by the shocks associated with the early stages of 
star formation, in M82 it is determined by the UV radia- 
tion prod uced by evolved stars and t he enhanced cosmic 
rays flux ijGarria-Burillo et alll2002ft . 

We are grateful to the technical staff in Pico de Veleta 
for their supporting during observations. This paper 
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TABLE 1 
Relative fractional abundances 



Molecule 




M82 




Orion Bar 


NGC 7023 


Horschcad 


TMC1 


L134N 




E. 


(0,0) 


W. 


IF 


PDR peak 


IR peak 


CP 




C 2 H (10 la cm" 3 ) 


58 


41 


30 


57(a) 




16 (e) 


7 . 2 U) 


15<«) 


c-C 3 H 2 /C 2 H 


0.02 


0.04 


0.02 


0.02 (b) 


. 03 (fc)(d) 


0.06 le) 


0.8 (9) 


0.33 (! " 


CH3C2H/C2H 


0.2 


0.4 


0.3 






<0.07 (e) 


1.4") 


0.07 (s) 


c-C 3 H 2 /HC 3 N 


> 6 


> 4 


> 3 


>1 W(<0 




19<=> 


. 4 (h) 


3(3), (i) 


CN/HCN 


4 


6 


4 


3(a) 


4-8^ 




0.07 (/ ' 


0.07 (,) 


CN/HC3N 


>120 


>58 


>80 


> 180< a > 






0.05 (/)(,l> 


0.9« 


HCO+/HOC+ 


35 


49 


50 


<166 <a) 


50-120 (6) 
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Fig. 1. — Observed spectra towards the positions East (E), Center a nd West (W) in M82. The 3 0m beam at 90 GHz around the three 
observed positions has been drawn in the interferometrid HCO image by Garci'a-Burillo ct al. (2002). We have fully sampled the HCO disk 
at 3mm. 
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Fig. 2. — Predictions for the abundances of various species derived from updated Le Bourlot et al.'s model. The calculations have been 
carried out for n=4xl0 5 cm -3 and Go=lxl0 4 in units of the Habing field. The cosmic ray flux is set to £=5xl0 -17 s _1 (galactic value) 
in left panels and f=4xl0 — 15 s — 1 (M82 value as derived by Suchkov et al. 1993) in right panels. We have shadowed the region of the plot 
in agreement with the observational results in M82. 



